Proposed fitting function for the critical Casimir force on 4 He film below the A point 
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We have extended the mean field calculation 1] of Zandi et al and have obtained an approximate 
mathematical expression for the Casimir scaling function #(y), which if extrapolated to the domain 
ty 2 > —y > 0, becomes surprisingly similar to that obtained experimentally |2| by Ganshin et al. 
The extrapolated i?(y) can be regarded as a proposed fitting function, that appears to agree better 
with the experiments [2J in D=3, than the exact mean field result of Zandi et al [lj. 
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1. Introduction 

It was proposed two decades ago by Nightingale and 
Indekeu that, the confinement of the critical fluctuations 
may give rise to a (classical) Casimir force [3J. Thereafter 
an experimental verification of the same by Mukhopad- 
hyay and Law [4j made the study of the classical Casimir 
force on critical films a hunting ground to the experi- 
mentalists and theoreticians. On this issue, the Casimir 
effects on different critical films have been the subject of 
a number of experimental [j42| and theoretical [l|, llCH16J 
works within the last few years. 

Few years ago Garcia and Chan [5| and Ganshin et al 
investigated the temperature (T) dependence of the 
critical Casimir force, and measured the Casimir force 
induced thinning of the liquid 4 He film near the bulk 
A point (T\ = 2.1768 K). They obtained a universal 
scaling function (#) of the Casimir force, and observed a 
dip minimum in the Casimir scaling function below the A 
point. This experiment challenges our understanding of 
the finite size effects on the films near their bulk critical 
points. 

For T > T\, a renormalization group calculation for 
the $ of the 4 He film was nicely presented by Krech and 
Dietric h I l71| . Although the calculation of Krech and Di- 
etrich [I7l | matches well with the experiment [5j, yet no 
satisfactory theory was given so far for T <T\. However, 
the mean field calculation of Zandi et al |l| for this re- 
gion, represents the basic feature of the classical Casimir 
force. 

For T < T\, Zandi et al [l| obtained an analytic expres- 
sion of "d in terms of the maximum (0o) of the supcrfluid 
order parameter (0). It is necessary to know 0o for plot- 
ting the mean field part of the Casimir force. Although 
the graphical solutions of the 0o are exact, yet the graph- 
ical solutions do not appear in the closed form. With a 
motivation of getting an approximate expression of 0o in 
the closed form, we extend the calculation of Zandi et al 

Since the maximum correlation length (£) of a critical 
film (of thickness L) is L/n [1|, the order parameter ap- 



§ars into the film for £ < ^, and it vanishes for £ > ^ 
. Although our approximate mathematical expression 
of $ for £ < — , matches very well with the exact mean 
field result of Zandi et al [ll], yet we extrapolate our ap- 
proximate expression of I? to the region £ > — where 
the order parameter actually vanishes. This extrapola- 
tion surprisingly fits the experimental data obtained by 
Ganshin et al [2j, and appears to agree better with the 
experiments 2], than the exact mean field result of Zandi 
et al 0. 



2. Exact mean field result 

Before going into the details of our work, let us briefly 
reproduce the calculation of Zandi et al [1]. For T <T\, 
we start from the Ginzburg-Landau free energy 



F = A 



dz 



+ -0 2 (z)+u(f> 4 (z)\dz, (I) 



where r = — 1/£ 2 , £ = £ \T/T\ — l\~ l/ , v is the correlation 
length exponent, L (~ 300A [2j) is thickness of the film 
along the z direction, A (~ 1 inch 2 [2|) is the area of the 
film along the x — y plane, u is a positive coupling con- 
stant, and the mean field order parameter (4>(z)) obeys 
the equation [l| 



(*) 



dz 2 



+ r<j){z) +4u0 3 (z) = 0. 



(2) 



Solutions of the Eqn.(2) with the Dirichlet boundary con- 
ditions (0(0) = (p(L) = 0) arc the Jacobian elliptic func- 
tions [18j]. For the ground sate, there is a single maximum 
of 4>{z) at z = L/2. With the consideration of Dirichlet 
boundary conditions (0(0) = 4>(L) — 0), and that of 
a single maximum of (j){z) at Z = L/2, we can write a 
transcendental equation for the maximum of the order 
parameter (0(4) = 0o) as [l| 



VT=W=i//2 = A" (r?/(l -»/)), 



(3) 



where y 
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is the complete elliptic integral of 
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3. An approximate mean field result in a closed form 

Although the graphical solutions do not give a mathe- 
matical expression of r\ in terms of y/— y/2, yet we need 
a mathematical expression of r\ for conveniently plotting 
the Eqn.(3). This purpose can approximately be solved 
by interpolating the asymptotic solutions near r\ = and 
77 — > |. It is easy to check from the Eqn.(3) that, the 
allowed range of 77 is < 77 < \. For 77 — > 0, - v / 3 y/2 goes 
to 7r/2, and we can expand the 77 parts of the Eqn.(3) up 
to the third order in 77 to get 



FIG. 1: The dots are the graphical solutions of the Eqn.(3). 
The dotted, dashed and continuous lines follow the Eqns.(5), 
(6) and (7) respectively. 
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first kind. Since, the integration limit and 4>{z), both are 
functions of L, the derivative of F (in the Eqn.(l)) with 
respect to L is nontrivial. The mean field force acting 
on the film, can be shown to be — |r = ^-h(ij) I r 

' 5L 2 V dz / I z—L 

[l(. Subtracting the bulk force (— 337 1 L . ) from — |^, 
we can obtain the mean field Casimir force, which can be 
scaled by L 4 /(^ AfcsT\) to obtain the Casimir scaling 
function in the dimensionless form as [lj] 



<%) 



?r 



4^ uk B T) 



(\- f](V=v/2) (l-»7(V z i//2))). 

(4) 



To plot the scaling function (in the Eqn.(4)) we must 
need to know rj{y/—y/2) from the Eqn.(3). Since 77 by 
the definition is positive, one can check from the Eqn.(3) 
that, a non zero positive r](y/—y/2) exists only for [1] 
—y > ir 2 . Although Zandi et al [if plotted the Eqn.(4) for 
00 > — y > 0, yet they did not clearly describe how they 
obtained rj(y/—y/2) for 00 > — y > ir 2 . However, one can 
easily check that rj(y/~y/2) be exactly solved from the 
Eqn.(3) by the standard graphical method. Within this 
method, we can plot the left and right hand sides of the 
Eqn.(3) with respect to 77 keeping a given yj— y/2 fixed. 
The intersecting point of the two plots is a solution (77) 
with respect to the given y/— y/2. In this way, we can get 

the set {77} for a given set {-^-j^}- We plot these sets of 
points in the Figure 1. Adjusting the depth (#(—7.39) = 
— 1.3±0.03 [2j) of the experimental minimum at y = — 7r 2 
with a suitable value of u, one can plot the Eqn.(4) (with 
the graphical solutions) to the domain 00 > — y > it 2 for 
obtaining the t)(y) almost similar to that of Zandi et al 
Q. For 7T 2 > -y > 0, Zandi et al (U plotted the Eqn.(4) 
by putting rj(^y/2) = 0. 



In the other asymptotic limit, i.e. for r\ — > i, ^J—y/2 
goes to infinity. In this limit, 77 can be expressed as 77 = 
|(1 - 5), where 5 -> 0. In the Eqn.(3), K(r)/(1 - 77)) 
can be written (by the definition) as an integral of the 
form L . p which tells that p = 

° v /p(2-p)(l-^(l-2p+p2)) 

has a logarithmic divergence, so that most of the integral 
would come from p — > 0. In this asymptotic limit, the 
Eqn.(3) can be approximated with the first order in S 



and p, as *-^- = f 
asymptotically gives 



dp 



ri(V=y/z) 



-coth z 



y/MS+p) 



V^y/z 



V2sinh 



for 



[7?]' which 



> 1. (6) 



Let us now find a smooth function for 77 as an in- 
terpolation for the whole range < 77 < |, in such a 
way, that, it fits to the two extreme ends of 77. Since 

4>{z) = */=£-tanh[zy/—r/2] is a solution of the Eqn.(2) 
for an unbounded situation, we can take a trial function 
from the asymptotic solutions (Eqn.(5) and (6)), that 0o 
for the bounded system considered by us, would be close 

y/(-y-^)/2^ 



to 



T^tanhf 



With this consideration, we 



can take the trial interpolation 



ViV^y/2) = -tanh 2 



for the whole range < 



VH/-O/2 



(7) 



?— < 00. In the Figure 1, 
we see that, the trial interpolation in the Eqn.(7), fits 
very well with the asymptotic solutions (f]{s/—y/2)) in 
the Eqn.(5) and (6). The trial interpolation also fits very 
well with the exact graphical solutions of the Eqn.(3). 
Hence, the approximate Eqn.(7) can be used for plotting 
the Eqn.(4). 
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FIG. 2: The continuous part (— y > tv 2 ) of the plot follows 
Eqn.(4), and the dotted part (tv 2 > — y > 0) is the extrapo- 
lation of the same equation. The depth of the experimental 
minimum [2] has been fitted for ... . \ — =- = 0.01237. A 
few number of experimental points are replotted from the 
Ref.Q] with the substitutions y = x/£ /v , v = 0.67016 @] 
and £ = 1.2 A [l], [ig|. The dashed line is redrawn like that 
obtained by Zandi et al [lj. 



4. Extrapolation of the mean field result 

Since the minimum of the experimental d is obtained 
at y = -7.39 ± 0.61 [E Q, and not at y = -ir 2 Q|, 
we extrapolate the Eqn.(4) (by sacrificing the fact that 
rj(y/—y/2) — 0) to the domain n 2 > — y > 0, with the 
motivations, whether the Eqn.(4) has a natural minimum 
near y = —7.39, and whether the experimental points be 
fitted with the extrapolation. 

Putting the Eqn.(7) in to the Eqn.(4), we extrapolate 
the Eqn.(4) to the domain 7r 2 > — y > for obtaining 
the Casimir scaling function (in the Figure 2) similar to 
that obtained experimentally in the Ref. [21 ■ Now, we see 
in the Figure 2 that for £ = 1.2A 0, [h|, Eqn.(4) has 
a natural minimum at y = —7.03, which is closer to the 
position (y = —7.39 ± 0.61 [l|, [2J) of the experimental 
minimum. The extrapolated part surprisingly fits the 
experimental data of Ganshin et al [2j , and it appears to 
fit better than the exact mean field result of Zandi et al 



:n. 



5. Conclusions 



The extrapolation of the mean field result to the do- 
main 7T 2 > — y > is questionable due to the fact that, rj 
is not being equal to zero, rather being negative in this 
domain. For this reason, the extrapolation does not form 
a basis of an improved theoretical approach. The extrap- 
olated part of the Casimir scaling function becomes sur- 
prisingly similar to that obtained experimentally 2] by 
Ganshin et al. Our extrapolated part can be regarded as 



a proposed fitting function, that appears (in the Figure 
2) to agree better with the experiments [2|in D=3, than 
the exact mean field result of Zandi et al [1] . 

It is to be mentioned that, the mean field result of 
Zandi et al was improved by themselves [1| only at y = 
by introducing the contribution of the critical fluctua- 
tions. The mean field result was also improved in the 
Ref. [20j by the confinement of the Goldstone modes and 
that of surface fluctuations. 
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